Parity nonconservation in heavy atoms: The radiative correction enhanced by the 

strong electric field of the nucleus 
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Parity nonconservation due to the nuclear weak charge is considered. We demonstrate that the 
radiative corrections to this effect due to the vacuum fluctuations of the characteristic size larger 
than the nuclear radius ro and smaller than the electron Compton wave-length Ac are enhanced 
because of the strong electric field of the nucleus. The parameter that allows one to classify the 
corrections is the large logarithm ln(Ac/ro). The vacuum polarization contribution is enhanced 
by the second power of the logarithm. Although the self-energy and the vertex corrections do not 
vanish, they contain only the first power of the logarithm. The value of the radiative correction is 
0.4% for Cs and 0.9% for Tl, Pb, and Bi. We discuss also how the correction affects the interpretation 
of the experimental data on parity nonconservation in atoms. 
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Atomic parity nonconservation (PNC) has now been 
measured in bismuth H , lead || , thallium || , and cesium 
H . Analysis of the data provides an important test of the 
Standard Electroweak model and imposes constraints on 
new physics beyond the model, see Ref. ||. The analysis 
is based on the atomic many-body calculations for Tl, Pb, 
and Bi 0] and for Cs [00] • Both the experimental and 
the theoretical accuracy is best for Cs. Therefore, this 
atom provides the most important information on the 
Standard model in the low energy sector. The analysis 
performed in Ref. 0] has indicated a deviation of the 
measured weak charge value from that predicted by the 
Standard model by 2.5 standard deviations a. 

In the many-body calculations the Coulomb in- 

teraction between electrons was taken into account, while 
the magnetic interaction was neglected. The contribu- 
tion of the magnetic (Breit) electron-electron interaction 
was calculated in the recent papers [p|pT|]. It proved to 
be much larger than a naive estimate, and it shifted the 
theoretical prediction for PNC in Cs. As a result, the 
deviation from the Standard model has been reduced. 
The calculations have already been used to get 

new restrictions on possible modifications of the Stan- 
dard model, see, e. g., Ref. jyj. The reason for the 
enhancement of the Breit correction has been explained 
in Ref. |0| . In the case of the Coulomb residual interac- 
tion the effect of the many-body polarization is maximum 
for the outer electronic subshell and quickly drops down 
inside the atom |3|-^| . The Breit interaction is more sin- 
gular at small distances than the Coulomb one. Hence, 
the polarization is maximum for the lowest subshell (Is 2 ) 
and quickly drops down towards the outer shells. The 
estimate of the relative effect of the magnetic polariza- 
tion gives Za 2 instead of naive a 2 , where Z is the nu- 
clear charge and a is the fine structure constant. To find 
the Breit correction there is no need to repeat the in- 
volved many-body calculations performed in Refs. || 
Indeed, the Breit correction comes from small distances, 



r ~ as I 'Z {as is the Bohr radius), while all the Coulomb 
polarization and correlation corrections come from large 
distances, r ~ ag. Therefore, it is sufficient to calculate 
the relative Breit correction to some PNC mixing matrix 
element (say 6s 1/2 — 6pi/2 mixing in Cs) in the simplest 
Hartree-Fock or RPA approximation. The relative Breit 
correction to the PNC effect with account of all many- 
body Coulomb polarization and correlation corrections is 
exactly the same. 

The Breit correction to PNC is just a part of the ef- 
fect. This part is related to the virtual excitations of 
the Is 2 subshell. Another contribution comes from the 
vacuum fluctuations, i.e. from the radiative corrections. 
Attempts to estimate this effect were made in Ref. pl| 
and gave very small values of the corrections. It has been 
pointed out recently [Q that the strong electric field of 
the nucleus enhances the radiative corrections, and they 
may be comparable with the Breit correction. Very re- 
cently this suggestion has been confirmed by the numer- 
ical calculation of the vacuum polarization correction in 
Cs 0. 

In the present paper we consider PNC in heavy atoms 
and calculate radiative corrections enhanced by the 
strong electric field of the nucleus. We calculate ana- 
lytically the leading term in the correction and we also 
estimate other terms. It turns out that in Cs, Tl, Pb, and 
BI the radiative correction compensates the Breit correc- 
tion calculated in Refs. [j9 10|. Thus, we return back to 
the result of the experimental data analysis made in Ref. 
P): deviation from the Standard model is 2.2 - 2.3 a. 

In the Standard model it is accepted to normalize the 
Weinberg angle at the W-boson mass My/. Atomic ex- 
periments correspond to a very low momentum trans- 
fer compared to My/. The renormalization from My/ to 
zero momentum transfer was performed in Refs. p^| , ^6[ . 
This renormalization is reduced to the logarithmically en- 
hanced single loop corrections cx a/ir that lead to what is 
usually called the radiative correction to the nuclear weak 
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charge. Account of this very important correction gives 
the nuclear weak charge Qw measured in on-mass-shell 
electron scattering at zero momentum transfer, where 
p 2 = m 2 and q = 0. However, atomic PNC corresponds 
to a different situation. The electron in the strong nu- 
clear electric field is off the mass-shell, p 2 ~ 1 jr 2 ^S> m 2 
and, besides, the typical momentum transfer is of the or- 
der of the inverse nuclear radius, q ~ 1/tq. In a precise 
calculation these effects must be taken into account. It is 
convenient to use Qw calculated in Refs. (ljjl^l as a ref- 
erence point. Then the renormalization procedure is the 
same as that in quantum electrodynamics and the correc- 
tion vanishes on mass shell at zero momentum transfer. 
In this approach it is clear that the correction we are talk- 
ing about is somewhat similar to the radiative correction 
to the hyperfine constant in a heavy atom [ fL7f . 

The wave function of the external electron is of the 
form 



u(r) 



F(r)n 
iG{r)Vt 



(1) 



where f2 and fl = — (er • n)Q are spherical spinors p8| . 
At small distances r -C ZaXc, where Ac is the electron 
Compton wave-length, the electron mass is small com- 
pared to the nuclear Coulomb potential, and the radial 
wave functions obey the equations 



d(rF) 

dr 
d{rG) 

dr 



n , , Za , „, 

-(rF) rG 

r r 



0, 



(2) 



K 



Za 



(rG) + (rF) = 0. 



For PNC effect we need to consider only Si/ 2 (« = — 1) 
and pi/2 (k = +1) electron states. Solution of Eqs. (0) 
reads 



F = Art-\ G = A- 



Za 



7 



(3) 



where 7 = yl — Z 2 a 2 and A is some constant dependent 
on the wave function behavior at large distances (r ~ as) 
fL9[ . In the leading approximation the PNC interaction 



related to the weak charge is due to Z-bozon exchange, 
see Fig. la. Calculation of the corresponding weak inter- 
action matrix element gives |1S[ 



< Pi/2\Hw\si/2 >o= Mo cx (F S G P - G s F p )\ r= 



(4) 



At ro — * this matrix element is divergent, Mq oc r 2 , 1 2 . 
As a result, the relativistic enhancement factor is R «3 
for Cs and R » 9 for Tl, Pb, and Bi U|. In the present 
paper we show that this divergence results in the double 
logarithmic enhancement of the radiative corrections. 

The first correction is shown in Fig. lb. It corresponds 
to a modification of the electron wave function because 
of the vacuum polarization. In the leading Za approx- 
imation the vacuum polarization results in the Uehling 
potential |20(| . At r <C Ac, this potential is of the form 
V(r) « 2Za 2 [ln(r/Ac)+C+5/6]/(37rr), where G ps 0.577 
is the Euler constant. Account of higher in Za cor- 
rections in the vacuum polarization leads to a modifi- 
cation of the constant: C — > C + 0.092Z 2 a 2 + see 
Rcf. plj] . However, this correction is small and can be 
neglected even for Za ~ 1. The potential V(r) mod- 
ifies the Coulomb interaction in Eqs. (||) —Za/r — > 
—Za/r + V(r). It is convenient to search for solu- 
tion of the modified Eqs. (g) in the following form 
T = F(l + F^), G = G(l + F^ 1 )), where F and G 
are given by (^) . The functions Fs]p and Gi^, satisfy the 
following equations 

7 + K dx 37T 

-J-H^-G<» + FW = -p X , (5) 
7 — k dx 3n 

where x = ln(A/r), and A = Acexp(— C — 5/6). The 
solution of these equations reads 



F« = — 
3tt 

Gi X) = ± 
3tt 



(Za) 2 2 + k 

7 7 Z 27^ 

(Za) 2 2 K (7 — K ) K 

7 7 2 27 2 



(6) 



where a is some constant. Using T and instead of F 
and G in Eq. (Q) we obtain the PNC matrix element in 
the form 



d e f 

FIG. 1. The leading contribution and the radiative correc- 
tions to the PNC effect. The double line is the exact electron 
Green's function in the Coulomb field of the nucleus, the cross 
denotes the nucleus, the zigzag and the dashed lines denote 
Z-boson and photon, respectively. 



< Pl/2 \H w \s 1/2 >=M (l + S), (7) 
where 5 due to the diagram Fig. lb is 

S b = ^ (fP + G«) + ^ (F« + <#>) . (8) 

To find the correction with the logarithmic accuracy 
there is no need to calculate a in Eqs. (0), it is enough 
to substitute the logarithmic terms from (^J) into fll) . An 
analysis that also includes a consideration of distances 
r ~ Ac gives 
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a [ — Za 
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2(Za) 2 

37T7 



[ln 2 {b\c/r Q )+B~ 



(9) 



where 6 = exp(l/(27) — C — 5/6), and B ~ 1 is some 
smooth function of independent of ro. A numerical 
calculation of 8b for Cs was performed recently in Ref. 
fuj . The result is in a good agreement with Eq. 
Comparison of Eq. (j^) with the result of Ref. |Q allows 
also to determine _B: We would like to emphasize 

that Eq. (H) does not assume that Za <C 1, it is valid for 
any Za < 1. Note that the (Za) 2 term in (^|) is larger 
than the Za one at Z > 10. 

We already pointed out that the weak charge calcu- 
lated in Refs. [p~5| , p~6| corresponds to zero momentum 
transfer. On the other hand, it is clear from Eq. (^) that 
the weak interaction matrix element is determined by 
the momentum transfer q ~ 1/Vo. The renormalization 
of the weak charge from q = to q ~ 1/Vo is described 
by diagrams c and d in Fig.l. A simple calculation gives 
the following correction SQw/Qw = 5 c d related to this 
renormalization 



4aZ 



(l-4sin 2 6W)ln(A c /ro) « -0.1%. (10) 



Where Qw ~ is the Weinberg angle, sin 6\y ~ 0.2230, 
see Ref. ||. Note that this correction is practically in- 
dependent of Z because Z/Qw ~ —Z/N w —0.7, where 
A is the number of neutrons. One can also obtain the 
correction ( fl0| ) using Eqs. (2a, b), and (3b) from Ref. 

Next we consider the contribution of the electron self- 
energy operator E. This operator being substituted to 
the Dirac equation, m — > to + E, leads to the Lamb 
shift of the energy level and to the modification of the 
electron wave function, see, e.g., Ref. fll8fl . As shown in 
Fig.le, this modification influences the matrix element 
of the weak interaction. The diagram Fig.le is not in- 
variant with respect to the gauge transformations of the 
electromagnetic field. However, the sum of the diagrams 
Fig.le and Fig. If (the vertex correction) is gauge invari- 
ant. It is convenient to represent the self energy operator 
as a series in powers of the Coulomb field of the nucleus, 
E = E + Ex + E 2 + see Fig.2. 
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FIG. 2. The electron self energy expanded in powers of the 
Coulomb field. The solid line is the free electron Green's func- 
tion, the cross denotes the nucleus, and the dashed line de- 
notes the photon. 

We need E(r, r'|e) at r ~ r' <C Ac, e ~ m. In this limit 
a calculation in the Feynman gauge with logarithmic ac- 
curacy gives 



E = p— ln(p /to ), 



Ei = 



Za 2 
47rr 



70 ln(p 2 /TO 2 ), (11) 



where p = p^^ffi, 7^ is the Dirac matrix, and p^ is the 
momentum operator. All the higher terms are not log- 
arithmically enhanced. Further calculation in Feynman 
gauge is rather involved because the diagram Fig. If is 
also logarithmically enhanced and there is a delicate can- 
cellation between a part of the E-contribution and the 
logarithmic part of Fig. If. To avoid all these complica- 
tions it is convenient to use the Landau gauge. In this 
gauge neither E„ (n = 0, 1,2...) nor the contribution of 
Fig. If contain the large logarithm ln(Ac/f*o)- Therefore, 
after the renormalization, taking zero momentum trans- 
fer as a reference point, the contribution of Fig. If is of 
the form Sf ~ Za 2 (l + fZa/n + ...), where we expect 
/ ~ 1. Since E does not contain logarithms, at r <C Ac 
we have Ew = (Za 2 /irr)T>u, where T> ~ 1 is some matrix 
dependent on Za, and u is the electron wave function. 
Substitution of Ew into the Dirac equation results in the 
x-independent terms of the order of ~ a /it in the right 
hand sides of Eqs. (^). As a result, the relative correction 
to the matrix element of the weak interaction due to the 
diagram Fig.le contains logarithmically enhanced Z 2 a 3 
terms. There is also a Za 2 contribution coming from the 
distances r ~ Ac- All in all the total contribution of 
diagrams Fig.le and Fig. If is of the form 



Za 1 



Za . 
a x +a 2 ln(Ac/-r ) 

7T 



(12) 



Our preliminary estimate gives di ~ 0.15, and we expect 
a 2 ~ 1. Therefore, the value of 5 e f for Cs is 8 e f ~ 0.1%. 
The calculation of the coefficient a 2 in Eq. ( |l2| ) is a very 
interesting and challenging problem. At this stage we 
can claim only that the contribution ( |l2] ) is much smaller 
than that of the Uehling potential (|9|) because it does not 
contain the big logarithm squared. 

There is one more radiative correction that has never 
been considered before. This contribution is due to the 
virtual excitation of the nuclear giant dipole resonance 
A* shown in Fig. 3. 




FIG. 3. Virtual excitation of the giant nuclear resonance. 



Our estimate gives 



-0.1Z 2/3 a 2 . 



(13) 



So it is completely negligible. 

Considering all the corrections, one has to remember 
about the contribution of the electron-electron weak in- 
teraction. Although this is not a radiative correction, we 
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still denote it S e - e . According to Ref. |2^] this correction 
is 



0.56 
ZR{Z) 



(1-4 sin 2 W ), 



(14) 



where R(Z) is the relativistic enhancement factor defined 
in Ref. [^9| . Values of R are presented after Eq. ([!]) . For 
Cs the <5 e _ e -correction is +0.04%, so it is also negligible. 

In summary, we have calculated radiative corrections 
to the parity nonconservation effects in atoms due to 
the strong electric field of the nucleus. The structure of 
the correction is somewhat similar to that for the radia- 
tive correction to the hyperfine interaction [ p7| . However 
there are two essential differences. 1) The corrections to 
the PNC effects are enhanced by the large logarithm de- 
pendent on the nuclear radius. 2) The contribution of the 
Uehling potential to the PNC radiative correction is dom- 
inating because it contains the second power of the log- 
arithm, see Eq. (^). The vertex and the self-energy cor- 
rections do not vanish, but they are not enhanced by the 
second power of the logarithm. The vertex, self-energy, 
nuclear excitations, and the electron-electron corrections, 
see Eqs. ©, (||), (|f), and ([[J), contribute at the level 
±0.1%. 

The correction (||) has been derived for the electron- 
nucleus weak interaction independent of the nuclear spin, 
i. e. for the weak charge Qw- However, the same formula 
is valid for the PNC effect dependent on the nuclear spin, 
i. e. for the anapole moment. This is obvious because the 
anapole moment interaction is also local, see Ref. JTsf ] . 

The radiative correction @ is +0.4% for Cs and +0.9% 
for Tl, Pb and Bi. The value of the Breit correction 
in Cs is —0.6%, see Refs. The Breit correction 

scales practically linearly with Z, see Ref. |0|. Therefore, 
for Tl, Pb, and Bi the value of the Breit correction is 
—0.9%. Thus, the radiative correction compensates the 
great part of the Breit correction in Cs, and practically 
exactly compensates it in Tl, Pb, and Bi. With account 
of both corrections, the difference between the Standard 
model prediction and the experimental result in Cs, see 
Ref. [§), stays at the level of 2.2a. 

Concluding we would like to say that both the radia- 
tive and the Breit corrections presently are known with 
high accuracy and this accuracy can be even further im- 
proved. The present limitation of the accuracy in the 
atomic calculations comes from many-body effects with 
pure Coulomb electron-electron interaction considered in 
Refs. An improvement of the accuracy of these cal- 

culations is a feasible but a very challenging problem. 
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